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Abstract 

For the first time, a full one-loop calculation of the process XX 27 is performed, 
where x is the lightest neutralino in the minimal supersymmetric extension of the Stan¬ 
dard Model. This process is of interest for dark matter detection, since it would give 
a sharp 7 ray line with Ej = m^. We improve upon and correct published formulas, 
and give cross sections for supersymmetric models with x masses between 30 GeV and 
several TeV. We find a new contribution, previously neglected, which enhances the 27 
rate for TeV higgsinos by up to an order of magnitude. As a byproduct, we obtain a new 
expression for the related process xX ^ 2 gluons, which on the other hand is generally 
smaller than previously calculated. 

There has been a recent claim that evidence for a 3.5 TeV higgsino annihilating into 
a 7 line may already exist from balloon emulsion and Air Cherenkov Telescope data. 
We comment on attractive features and problems with this interpretation. 


* lbe@physto.se 
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1 Introduction 


Most extensions of the Standard Model of particle physics contain supersymmetry in one or 
another form. In particular, if gravity is unified with the other fundamental forces, it seems 
that supersymmetry is a necessary ingredient. The known particles of nature seem to be 
unstable, with very few exceptions. Of course, the stable ones have a specific significance 
since they have survived since the times immediately following the Big Bang, and today form 
the visible structures (clusters, galaxies, stars, planets...) in the Universe. 

If there exist stable particles in the supersymmetric sector, those should also have survived, 
and could give a contribution to the total mass density of the Universe in a way which is 
calculable once the supersymmetric model is specified. Models which could contain stable 
particles are usually endowed with a symmetry, i?-parity, which guarantees that the lightest 
supersymmetric particle is stable (for a review, see [Q). If this particle is electrically neutral, 
as is the case for the neutralino (a linear combination of the superpartners of the photon, the 
Z boson and two neutral CP-even Higgs bosons), its weak interaction strength provides a 
relic density which can be within an order of magnitude of the critical density of the Universe, 

Peril = 1-9 • g/cm^, (1) 

where h is the Hubble parameter in units of 100 kms“^Mpc“^. Present measurements H, ^ 
indicate that h is in the range 0.5 — 0.7. It is customary when discussing energy density to 
normalize to the critical density through the dimensionless quantity H: 

H=-^. (2) 

Perit 

With the present value of h, the critical density H = 1 corresponds to = 0.25 — 0.5. 

From dynamical estimates of the total mass density on cluster and supercluster scales, it is 
found that Qtot > 0.4 [Q. An analysis based on supernovas of type la as standard candles 
Q is not inconsistent with H = 1, a value which has a strong theoretical preference since 
it is predicted in generic models of inflation. The contribution Hf, from ordinary matter 
(mainly baryons) is limited by nucleosynthesis arguments to be less than around 0.1 |^. 
Also, structure formation from the primordial fluctuations seen in the microwave background 
seems to be problematic unless there is a considerable amount of dark matter 0 . 

The “missing mass” problem found in the analysis of galactic rotation curves [|| may thus 
have as its solution the existence of a relic population of neutral massive particles residing in 
the galactic halos. We will assume that these particles are neutralinos x = Xi^ the lightest 
supersymmetric particle in the minimal supersymmetric extension of the Standard Model. 


2 Detection methods 

Several methods have been discussed to detect neutralinos (for a comprehensive list of ref¬ 
erence for all methods mentioned here, see Q). One possibility is to register scattering of 
neutralinos directly in sensitive terrestrial detectors, well shielded from radioactive and other 
backgrounds. After several years of development, some of these experiments now have a 
sensitivity large enough to rule out some supersymmetric models which otherwise satisfy all 
accelerator and other constraints |^, |l^ . 
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Other, more indirect, methods utilize the fact that neutralinos in the halo (or those accumu¬ 
lated gravitationally in massive bodies like the Earth and the Sun) may annihilate, resulting 
in ordinary particles which may be detected. For example, antiprotons and positrons which 
are not produced in large quantities by cosmic rays, would be produced in annihilations in 
the halo and would be continuously injected and “stored” during around 10 ^ years which is 
the containment time of cosmic rays in the galaxy. Although the fluxes could be of measur¬ 
able size the problem is that they are generally not much larger than the cosmic-ray induced 
fluxes, which necessitates a difficult model-dependent background subtraction. Also, since 
the spatial and energy distributions of these antiparticles are rather featureless, the lack of 
definitive signature is a problem. This is also true for annihilation into particles which give 
secondary gamma rays. Although this could give a broad feature in the spectrum up to the 
neutralino mass, it is in most supersymmetric models too weak to be convincing || 0 1 . 


A much better signature is given by neutralino annihilation into neutrinos near the centre of 
the Earth or the Sun. Since the ordinary solar neutrinos only have at most MeV energies, 
a multi-GeV neutrino signal from the Sun (or the Earth) would be an unmistakeable signal. 
If the neutrino detector has an angular resolution better than a few degrees, the mass of the 
annihilating particles can be determined with reasonable accuracy [Q. Also this method 
has already given useful constraints on supersymmetric models of dark matter ||^, and has 
a large future potential as large new neutrino telescopes are being built 0 - 


Finally, an excellent signature is provided by the annihilation of neutralinos in the halo 
into two-body final states containing a photon, such as 77 . This process was first proposed 
by Bergstrom and Snellman ||^, who also_^rst estimated the related two-gluon process (for 

ITT “ 


photinos). Other early references include |16 

10-3 


|18| . Since the annihilating neutralinos move 


with galactic velocities, r^/c ~ IQ-^, the result is a gamma ray line of phenomenal peak-to- 
width ratio. This signature has no conceivable background from known astrophysical sources. 
To optimally utilize this sharp signal peak, detectors with energy resolution at the percent 
level and below should ideally be used. In the region 30 - 100 GeV, where “light” neutralinos 
(mostly rather pure binos) are likely to reside, a detector has to be placed in orbit (above 
the atmosphere) not to be flooded with cosmic-ray produced atmospheric gammas. This led 
to a (temporary?) impasse for this idea, since photon detectors with high energy resolution 
are generally very heavy and therefore costly to launch. 


It was subsequently realised, in work by Urban et al , that a certain type of ground-based 
detector (Air Gherenkov Telescopes, ACT) could be used for searching for gamma lines from 
very heavy neutralinos (several hundred GeV to the TeV region). Here the limited energy 
resolution (10 %, at best) is partly offset by the very large collection areas possible ( 10 ^ to 10 ® 
m^). The rates for TeV neutralinos (generally higgsinos, since binos of such masses would 
overdose the Universe) were calculated in |^. In this paper, the process XX ^ ^7 was 
also proposed as an important channel, and a lower bound based on unitarity for high-mass 
higgsino annihilation into 77 and Z 7 was computed. In [pT| , a more complete calculation of 
the 77 process was presented, and the phenomenology was also worked out in (although 
in none of these TeV neutralinos were considered). 


Since there is a new generation of Air Cherenkov telescopes being planned |^, and also a 
new ambitious satellite project we have considered it timely to revisit the 77 process, 
performing for the first time a full one-loop calculation improving on the estimates of previous 
work. As a byproduct, we also recompute the two-gluon process, correcting a numerically 
significant error in . We give analytical expressions for all the amplitudes and give results 
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of scans in supersymmetric parameter space for several sets of models. 

During the final stages of this work, we have come upon a recent paper where it is 
claimed that evidence for a TeV higgsino annihilating into a 7 line may already exist from 
balloon emulsion and ACT data. We will comment on attractive features and problems with 
this interpretation. 

It is important to note that all the different methods to detect supersymmetric dark matter 
are complementary in many respects. Each process has certain regions in parameter space 
where it dominates. In fact, for TeV higgsinos, the 7 line process may probe supersymmetric 
models out of reach for accelerators in the forseeable future (for other possible methods to 
detect heavy higgsinos, see p^). 


3 Cross section 


We compute the full expression for the annihilation cross section of the process 
A? + X? ^ 7 + 7 


( 3 ) 


in the limit of vanishing relative velocity of the neutralino pair. This result can be applied to 
neutralinos annihilating in the galactic halo, where it is estimated that massive dark matter 
particles move with a velocity of the order u/c ~ 10“^. The outgoing photons will then be 
nearly monochromatic, with an energy equal to the mass of Xi- Furthermore, for very low 
relative velocity, the neutralino pair must be in a S wave state with pseudoscalar quantum 
numbers; we can therefore project out of the amplitude the ^Sq state, using the projector 


Ops 



( 4 ) 


where is the mass of Xi and p = 0) its momentum. Analyzing the Dirac structure of 

the Feynman diagrams we can draw for this process, we can exclude a priori those diagrams 
which cannot generate a pseudoscalar state. 

As suggested in Ref. it is convenient to compute the W boson and unphysical Higgs 
loop diagrams choosing the non linear gauge defined in Ref. |^, which is a slight variant 
of the usual linear R-gauge (or ’t Hooft gauge). In this particular gauge the 7 

vertex vanishes and some of the box diagrams which, in the linear R-gauge, may give a 
contribution to the amplitude are not present. In addition the 7 vertex and the W 

propagator assume simpler forms. We find that the Feynman diagrams we have to include 
in the computation, at one loop level, are those shown in Fig.|^. 

There is another property that introduces a further main simplification. In any process 
of annihilation of particles at rest, the four point loop integrals one obtain from each box 
diagram applying the Feynman rules, can always be rewritten as a linear combination of three- 
point integrals. This is possible because, for each integral, one can find a linear combination 
of the four factors in the denominator which is independent of the momentum flowing in 
the loop. It follows that we do not have any four-point integral in the computation of the 
amplitude. 
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The amplitude of the annihilation process can be factorized in the form 


6 ~ 

= n /?; , £(£i>e 2 ,fci,fc 2 ) A (5) 

zyz TT^ 

where ei, £2 and ki, k 2 are respectively the polarization tensors and the momenta of the two 
outgoing photons. The cross section is then given, as a function of A, by the formula 


va2-y 


a'^Ml 

IbTT^ 



( 6 ) 


We write the total amplitude as the sum of the contributions obtained from the four different 
classes of diagrams: 


A — Aj^ Aff+ + Aw + Ag ■ 

The single contributions are listed below. In each of them we separate real and imaginary 
parts. To facilitate a comparison, we follow as closely as possible the notation of pl| . 

1) Contribution of the fermion-sfermion loop diagrams (la-ld in Fig.|^. The sum over / 
includes the quarks and the charged leptons, the sum over / the corresponding sfermion in 
both chiral states. 


Re A 


ff 


Ml 


ff 


\ b S r f '\/ Qj b D r f 1 S r f 

9 + 9 
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(7) 
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Gzf 
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■ log 


f 1 + ~ 
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( 8 ) 


where mz and m^o are the masses of the Z boson and of the Higgs pseudoscalar H^, and 
we have defined: 


a = 


M} 


b 


2 

mi 

Ml 


Q - -L (M nb* nb*] 

^ff 2 r//i %/i ^^//i %i/ 


Dff = - 

ff 2 


(^//i ^//i + ^//i ^//i) 
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Gzf = 


sT] 

COS 9w 


(sill flfii) GijOf — 9 h°ii) {dHlff 9 hI 


// 


The index 1 is referred to 6/ fti® charge of the fermion in unit of the electron charge 
(—e), Cf is the color factor equal to 3 for quarks and to 1 for leptons, is the third component 
of the weak isospin equal to +1/2 for f=u,c,t and to -1/2 for f=e,/r,r,d,s,b. 

2) Contribution of the chargino-Higgs boson loop diagrams (2a-2d in Fig.||). The sum over 
xt includes the two chargino eigenstates. 


Re — 


T — 


1 b S^H + D^h j , , N ,1 j , 

^ - Ji (a, 5) + - h (a, b) 


2 1 + a — 5 

b Sy-H + Vo,b 1 b Sy-H 

1 + a — b 21 — 6 

M+M^Gho^ 


Am 


I 4 ( 4 M 2 -m^o) 



(9) 


ImAH+ = 0 

where we have defined: 

f 2 M ^4 


a = 


Mi 


'‘H+ 


b = 


^H+ 


( 10 ) 


^xH — 2 (Sff+li 9H+li + 9H+li Sff+li) ^ (Sff+li flff+li + 9H+li Sff+li) 


Gzx — {gzii 9zii) (szii sill) Gh°x ~ (^5“ii 9h°ii) {shIu 9H°ii) 

and the index i is referred to xf . 

3) Contribution of the chargino-W boson loop diagrams (3a-3c in Fig.^. 
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( 12 ) 


where in this case: 


M\ 

^ 


b = 






^xw — 2 {dwii 9w*ii + 9wii 9w*u) ^xw — 2 {9wu 9wii + 9wii 9w*ii) 
4) Contribution of the chargino-unphysical Higgs loop diagrams (4a-4b in Fig.|^). 


Re A-g = ^ ^ 


M 2 


1 Sgx + x/a Dg^ 

2 l-a -6 


h (a,l) 


+ 


1 Sgx Sgx + x/a Dgx \ r i' ^ ^ 


2 1-6 


1 — a — 6 


l 2 {a,b)-- 


2 1-6 


/a (a, 6 ) 


(13) 


IrnAa = 0 


(14) 


where a and b are the same as those defined for the chargino-W boson contribution and 


Sgx = \ {ghu 9 gu + 9Gii 9 gu) ^Gx = ^ (sgu 9gu + 9gu 9gu) ■ 


The coupling constants for left and right chiral states, and g^, are written in the notation 
adopted in the PhD Thesis of Edsjo all of them are defined therein, with the exception 
of the unphysical Higgs couplings, which we rewrite from Ref. in the conventions of pC|]: 


L 

9Gli 


9Gli 


= — gsin/? 


= +g cos (3 


Zl^Va 


-{Z*^ + Zl,Unew)V,l 


ZisUii — \ -(Z 12 + Zii tan9iv)l7i2 


£i 


(15) 

(16) 


The functions /i(a, 6 ), 12 ( 0 , 6 ) and 13 ( 0 , 6 ), which arise from the loop integrations, are 
specified in Appendix A. R (a, b) is the well known three point function that appears in 
triangle diagrams; it is an analytic function of a/b. I 2 {a,b) and I 3 {a,b) are given both in 
form of integrals over a Feynman parameter x, and may be expressed in terms of dilogarithms. 
It is numerically very easy to compute them even in the integral form as, for any physically 
interesting value of the parameters a and b, the integrands are smooth functions of x. In the 
real parts there are apparently poles for a — b = 1 , 6 = 1 , a + 6 = 1 : we have verified that, in 
each formula, a cancellation among the different contributions takes place for each of these 
spurious poles and therefore we do not have to deal with any numerical problem. 
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The results for the imaginary parts have been checked applying the Cutkosky cutting rules 
to the corresponding diagrams. Eq. ( 0 ) reproduces equations (2.5) in Ref. and (2.5) 
in Ref. except for one small, but numerically significant, error in both of these.f] As 
we will briefly sketch in Section ^ this discrepancy may have led, for some SUSY models, 
to an overestimate of the cross section of the annihilation of neutralinos into two gluons. 
In Eq. hxing the overall sign through Eq. (||), we do not agree with the sign of the 

result quoted in Ref. This sign difference is numerically important in cases where the 
amplitudes are of the same order of magnitude. 


4 Results and Discussion 

We have calculated the yy ^ 27 rates for a broad selection of supersymmetric models. 
All experimental constraints on the parameters have been checked, and only non-excluded 
models have been treated. (For a detailed explanation of the method of computation and 
the parametrization of models, see [@.) 

The calculation of the relic density Q,h? has been done using the full machinery of exact 
treatment of resonances and thresholds |Q, loop corrections to masses and vertices 
and coannihilations |Q. We are only interested in models which have a chance to at least 
explain the galactic halos, so we keep in all the graphs presented here only models which 
satisfy 0.025 < VLh? < 1. 

The detection rate of galactic gamma rays from neutralino annihilations depends on two 
different parts. One concerns the particle physics calculation of the 27 cross section for a 
given set of supersymmetric model parameters, and also the calculation of the relic abundance 
of neutralinos in standard Big Bang cosmology. The second part concerns the distribution of 
neutralinos in the galactic dark matter halo, which is a difficult problem of galactic structure 
formation. Unfortunately, all detection predictions are very sensitive to the poorly known 
dark matter mass distribution in the halo, since the 7 line flux from a volume element dV of 
density is given by 




And? 


dV, 


(17) 


where d is the distance to the volume element and the factor of 2 comes because there are 
two 7 s per annihilation. To obtain the full rate, one has to integrate this expression over all 
the halo (or, for a particular direction, over the line-of-sight). 

We start by discussing the halo-model independent quantities, which can be used for any 
particular halo dark matter distributions. Later we discuss the expected rates for some 
examples of halo models. 

In Fig. ^ we show the prediction for the annihilation rate into two photons, va 2 'y = va{x^ + 
Xi ^ 7 + 7) for ^ broad scan in supersymmetric parameter space (the parameter range is 
identical to that of the “generous scan” of Table 2.2 in [Q |^). As can be seen, for the 

*In the formula (2.5) of |25|, we think that the argument of the denominator of the second logarithm 
should be b — a(l — x^) instead of 6 + a(l — x^). This is plausible on general grounds, since it is precisely the 
zero of this expression in the integrand v^ich gives an imaginary part when > mf. It seems this error 
has propagated to the computer code of |l|. 
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low-mass region the rate varies over many orders of magnitude. In the most favourable cases, 
it is of the order of a few times 10 “^® cm^s“^, whereas it can be as low as 10 “® of that value 
for other models. This wide spread of predictions is typical of supersymmetry scans also for 
direct detection and neutrino detection rates. Low 27 rates are obtained for binos and mixed 
neutralinos, whereas higgsinos give a larger rate. This is particularly clear at the high mass 
end, where the requirement forces most of the allowed models to have a large higgsino 
fraction, and due to the W diagrams the 27 rate is exceptionally large for TeV higgsinos ||2^ . 
In this scan, the squark masses were generally quite large, which explains why also at low 
masses the models with the highest 7 rates are mainly higgsino-like. 

In Figs.^ (a) and (b), we have made more extensive scans for the low- and high-mass regions 
respectively. In the low-mass case, Fig.^ (a), we have allowed a larger fraction of models 
with low squark mass parameter (down to 150 GeV), since the bino rate into 27 is dominated 
by box diagrams of type (la) and (lb). As can be seen, some models give rates of the order 
of 10“^® cm®s“^. Now, because of the lower squark masses, models which give the highest 
rates have a large gaugino component. 

In the high mass region, we have focussed on the mass range between 3 and 4 TeV because of 
the preliminary indications in | 2 ^ of a possible signal. As can be seen, over that mass range 
there are many models that fulhl all accelerator and other constraints, and which all give 
rather high annihilation rates. The rates are in fact much higher than the lower limit, based 
on unitarity, given in (which we reproduce). Investigating the various contributions 
shown in Fig. ||, we have found that the reason for the large value of the (real part of) 
amplitude is diagram 3c, which gives a contribution proportional to 


1 _ -^1 

-L 71 f‘2 


Jr 




In the pure higgsino limit, the lightest higgsino Xi ^^nd the lightest chargino xt almost 
degenerate, so when the neutralino is much more massive than the W boson this t-exchange 
graph becomes dominant. The approximate treatment in |2(| only works if M^o / is 
much smaller than unity. 

In supersymmetric models, the mass difference between xf and xt goes down with increasing 
higgsino fraction of the neutralino. Thus, the 27 rate should go up with decreasing gaugino 
fraction Zg. This is clearly visible in Fig.|^ (a). In fact, we can approximate the pure higgsino 
limit Zg ^ 0 by keeping only the diagrams 3a-3c and inserting Sw = Dw = <?2,/4 {gw is the 
usual weak gauge coupling constant). Doing this, we find the roughly constant value 

?;cr 27 (pure higgsino) « 1.2 • 10 “^® cm® (18) 

over the whole range above 1 TeV.[| 

In Fig. H (b) we show the ratio 

F^ = 2 (19) 

_ vcrww _ 

^Although we have not yet computed the full amplitude for XX (L- Bergstrom and P. Ullio, in 

preparation), we have checked in the same limit that that process also tends to a constant value of gv of the 
same order of magnitude. 
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for our set of sampled supersymmetric models, bracketed by the unitarity lower bound of 
p0| and the pure higgsino rate found here. Notice the remarkable fact that a pure higgsino 
of 3 - 4 TeV mass can have the 27 channel as one of the dominant annihilation modes.| 

It is an interesting coincidence that a higgsino-like neutralino in the 3-4 TeV range gives a 
value of the relic abundance which is just right to give closure density, ^ 0.5. To arrive 
at this conclusion, it is essential to include in the relic abundance calculation coannihilation 
on the near-degenerate states, something which has only very recently been done by Edsjo 
and Gondolo (see Fig. 3 of |3^). 


5 Halo models 


In the most conservative model, the dark matter halo is described as an isothermal sphere 
with a density profile of the form: 


P{r) = Po 


( a^ + R^ \ 
\a'^ + ) 


( 20 ) 


where r is the distance from the galactic centre, po — 0.3 GeV/cm^ is the local halo density, 
a the core radius and R ~ 8.5 Kpc is our galactocentric distance. The gamma ray flux in a 
given direction is ^ : 




Va2'y 


J qJ: sight 

4.10-10-^ 


p^(0 diii^) 


10 GeV 


IQ-^o cm^s-i 





cm ^s isr 1 


( 21 ) 


where ip is the angle between the direction of the galactic centre and that of observation. 
The function J, which is obtained performing analytically the integral over the line of sight, 
is decreasing in the interval [0,7r] and its slope and maximum value increase as the parameter 
R / a increases. A typical estimate of J in the direction of the galactic centre is between 
J(0,l/3) - 0.4 and J(0, 2) - 2. 

Several models of the galactic halo contain the hypothesis of an enhancement of the dark 
matter density in correspondence with the galactic centre, where there seems to be a massive 
nucleus. The isothermal model of Ipser and Sikivie [ ^ , in which is predicted a large increase 
of the density in a region of 150 pc around the galactic centre, has been considered in Ref. Q . 
The estimated flux in the direction of the galactic centre is of the order of 10^ where 
<I>-y = $.y(0) is the flux that would result in the conservative halo model choosing J = 1, and 
again the value is to a large extent dependent on the parameters in the model. 

In the Berezinsky et al.|^ model, the dark matter halo density has the following profile: 

/ D \ 1-8 

p{r)=pu(^-j (22) 

^To be fully consistent, this mode should then also be included in the determination of the relic abundance. 
This necessitates a calculation of the momentum-dependence of the cross section, which has not been done. 
The effect is not expected to be dramatic. 
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which is assumed to be valid within 0.01 pc from the center. In this case the gamma ray 
source would appear pointlike with 99% of the flux emitted concentrated in a central region 
defined by a solid angle of less than 10“® sr. The flux in the direction of the galactic centre 
is highly enhanced, up to lO”^ 

An attractive hypothesis is that strings and texture fluctuations may have generated per¬ 
turbations in the dark matter halo providing a clumpy density distribution [ ^ . Also, in 
standard CDM scenarios, smaller “clumps” of dark matter may merge to larger clumps like 
galactic halos, and some of these objects may survive tidal forces. A significant increase in 
the gamma ray flux can be obtained if typical 10® Mq clumps are present. It is very difficult 
to detect such dark matter concentrations, nevertheless some of these may already have been 
identified, in the form of dwarf spheriodals residing in the Milky Way halo. As recently as 
in 1994, the Sagittarius dwarf was discovered, a dwarf spheroidal that is at a distance of 
about 15 kpc from the centre of the Milky Way, well inside the galactic halo potential, and 
about 23 kpc from the Sun. It is dark matter dominated and has a density profile resem¬ 
bling a Heaviside function with a value of 1.5 GeV cm“® on a length scale of 1 Kpc. A nice 
confirmation of the hypothesis of a line signal from dark matter annihilation would be the 
detection of an increased flux in precisely the direction of the Sagittarius dwarf. The flux in 
this direction would be larger by a factor of around 2 than the flux from the galactic centre. 

Taking the uncertainty of halo models into account, we can write the detected gamma line 
flux at Earth as 




(23) 


where we have seen that k can vary over a large range, from around 1 to at least 10®. 

In Fig. H we show the normalized flux 4)..^ for the same models as those in Fig.^. 

It is clear that the flux reported in |2^, ~ 10“® cm“^s“®sr“®, cannot be explained by any 
of the scanned supersymmetric models in the conservative halo model. However, taking the 
pure higgsino rate ( [l8| ) in the model of Berezinsky et al, such a high flux can be obtained 
towards the centre of the galaxy. However, the (weaker) indications of a high flux in other 
directions claimed in are not explicable in models where the density is only enhanced near 
the galactic center. A clumpy distribution with a large number of clumps rather smoothly 
distributed in the halo would be needed to explain the high flux of |26 . 


6 The two gluon channel 

As a final application of our loop calculation, we consider the annihilation channel XX 
2 gluons. Here only the diagrams la-ld contribute, with of course no contribution from 
leptons. To get the cross section, one just substitutes the electric charge Cj with 1 in the 
expressions (^ and (^) above, and the colour sum and average is performed by the simple 
recipe [|^ > 2a^ in (^). In Fig.^ we show the value of va 2 g for the scan of low-mass 

neutralino models using our results compared to those obtained with previous formulas . 
As can be seen, our prediction is for most of the models below the previous ones, which could 
make the indirect method of antiproton detection more difficult than previously thought. 
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7 Conclusions 


We have reinvestigated the process xx —> 27 in considerable detail. Our new results supersede 
and correct previous partial results in the literature. A new contribution is identified, which 
has its main effect for heavy higgsino-like neutralinos. A TeV-scale pure higgsino may in fact 
have the 27 final state as one of its main annihilation modes, in addition to giving a relic 
density close to the critical density of the Universe. 

A clumpy halo is needed, however, to explain the high rates indicated in if that analysis 
would survive the addition of further observational data. Irrespective of the fate of that 
3.5 TeV higgsino the gamma line signature is a promising one for neutralino detection and 
should be kept in mind when designing new gamma ray detectors. For a smooth halo, the 
ideal gamma ray telescope should have large collection area ( 10 “^ m^ or larger) and good 
energy resolution (down to percent level) but need not have good angular resolution. A large 
angular acceptance would be much more important. 

As a byproduct, we have also recalculated the xx ^ ‘^9 process, where we agree with pub¬ 
lished results apart from a small, but numerically significant, difference. 
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Appendix A 


We define below the functions h (a, 6 ), I 2 (a, b) and I 3 (a, b) which appear in the formulas for 
the real part of the amplitude. They arise from the three point loop integrals. 


h {a, b) 



Aax"^ — iax + b 


(A.l) 


h {a,b) 



—a a;^ -I- (a -I- 5 — 1 ) a; -I-1 


-I- (—a -I- 5 — 1) a: -I- 1 


(A.2) 


h {a, b) 



—a x'^ + {a + 1 — b) X + b 
ax^ + {—a + 1 — b)x + b 


(A.3) 


Ii (a, b) is the well known function that appears in the calculation of the real part of a triangle 
diagram. It is possible to perform the integral over the Feynman parameter x, obtaining the 
analytical form: 


II 











-2 f 




arctan 


1 


VE/a^ 


Ii (a,b) = < 


log 


— b/a 
1 —y/l—b/a 


if b > a 


if b< 


(A.4) 


I 2 (a,b) and I 3 (a,b) can be expressed in terms of dilogarithms. Li 2 {x) is defined according 
to Lewin iHl. We find: 


I 2 {a, b) = — Li 2 

+ Li2 


h (a, b) = - Li 2 
+ Li2 



(A.5) 


(A. 6 ) 


where we have defined: 


Ai = (a + 6 —l)^+4a, A 2 = (6 —a—1)^ — 4 a 


Let us consider first equation |A.5| . The second member is explicitly real. Va, b > 0, we find 
that Ai > 0 and that the arguments of the first two dilogarithms (1 — a — 6 ± \/A7)/2 are 
real and < 1: the corresponding dilogarithms are therefore real. If A 2 > 0, Va, 6 > 0, the 
argument of the third and forth dilogarithm (1 + a — 6 ± V^ 2)/2 are again real and < 1. If 
A 2 < 0, they are complex conjugated and we can use the identity 

Im {Li2{z) + Li2{z*)) = 0 

to show that there is no imaginary part. In the latter case it is easier to compute numerically 
I 2 ia,b) introducing LiMr,9), the real part of the dilogarithm of a complex argument 
We rewrite equation |A.5| as: 


I 2 {a,b) 


1 — a — b — 


— Li2 


+ 2 Li 2 (c(a), 9{a, b)) 




— Li2 


1 — a — b + VAi \ 
2 


•where r{a) = ^/a and cos 9{a, b) = {1 — a — b)/2 ^/a. 


(A.7) 
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It is not necessary to implement numerically equation A.6. Exploiting Landen’s functional 
equation Q: 


Li2{z) + Li2 




it is possible to express J3 (a, b) in terms of I2 {a, h). We have found the following identity: 
I3 (a, b) = -I2 (a, b) + log2 -—- - log^ -—- (A.8) 


2Vb J V 2Vb 

For A2 < 0, the second logarithm squared has a complex argument and takes the form: 


- log2 ^ - log^ 


2 f\z\ exp{iip) 


= ^ 


where (f = arctan((l + b — a)/\/—/S.2) ■ 
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Figure 1: Feynman diagrams contributing, at one loop level, to neutralino annihilation into 
two photons. Diagrams with exchanged initial and final states are not shown. 
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Neutralino mass (GeV) 


Figure 2: Annihilation rate of neutralinos into 2 7 for a broad scan in supersymmetric 
parameter space (the parameter range is identical to that of the “generous scan” of Table 
2.2 in [34,30]). 
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Neutralino mass (GeV) Neutralino mass (TeV) 


(a) (b) 

Figure 3: Annihilation rate of neutralinos into 2 7 for more extensive scans in the low- and 
high-mass regions. 



Gaugino fraction Neutralino mass (TeV) 

(a) (b) 

Figure 4: In Fig. (a) ratio 2 / vaww as a function of the gaugino fraction Zg] in Fig. (b) 

ratio 2 u(T 27 / vcrww for the sampled supersymmetric models of Fig. (a), the unitarity lower 
bound of [ 20 ] (dashed line) and the pure higgsino rate obtained in this work (solid line). 
Adding loop corrections to neutralino masses [33] (not done in these plots) one may find 
points closer to the solid line. 
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